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We consider the low-energy effective field theory describing the infrared dynamics of non- 
dissipative fiuids. We extend previous work to accommodate conserved charges, and we clarify 
the matching between field theory variables and thermodynamical ones. We discuss the systematics 
of the derivative expansion, for which field theory offers a conceptually clear and technically neat 
scheme. As an example, we compute the correction to the sound-wave dispersion relation com- 
ing from a sample second-order term. This formalism forms the basis for a study of anomalies in 
hydrodynamics via effective field theory, which is initiated in a companion paper. 



I. INTRODUCTION 

Low energy effective field theory is an extremely pow- 
erful tool to describe the dynamics of various experimen- 
tally accessible physical systems and to parametrize our 
ignorance of short distance physics. Techniques based on 
effective field theory are especially successful when sym- 
metries determine the low energy particle content and 
interactions. A classic example of such a situation is 
provided by the chiral Lagrangian, describing Goldstone 
bosons such as pions. As a consequence of the symme- 
tries, interaction between Goldstone bosons are weak at 
energies small compared to the symmetry breaking scale, 
so that the derivative expansion is a natural perturbative 
scheme to describe their dynamics. 

A rich class of phenomena where the low energy de- 
grees of freedom and their interactions are also fixed by 
symmetries to a large extent, is provided by hydrody- 
namics. Like the chiral Lagrangian, hydrodynamics is 
naturally organized as a derivative expansion, with the 
mean free time and the mean free path playing a role 
similar to the symmetry breaking scale for the chiral La- 
grangian. The similarity with the pion chiral Lagrangian 
goes even further. Indeed, hydrodynamical degrees of 
freedom are actually Goldstone modes, either of space 
translations spontaneously broken by the presence of the 
medium (phonons) [l| , or of global conserved charges car- 
ried by the fiuid [3|- 

Nevertheless, the traditional approach to hydrodynam- 
ics is quite different from conventional effective field the- 
ory. One starts with a set of conservation law for a set of 
"composite" objects — the energy-momentum tensor and 
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the conserved currents. The derivative expansion enters 
at the level of the "constitutive relations," which express 
the composite objects in terms of more elementary fluid 
quantities — the fiuid velocity field and local thermody- 
namical variables. As discussed in more detail below, 
this prescription becomes somewhat ambiguous as one 
goes to higher orders in the derivative expansions. 

The main goal of this paper is to demonstrate that, at 
least for a fiuid without dissipation, hydrodynamics can 
be recast into the traditional effective field theory lan- 
guage. Namely, just like for the pion chiral Lagrangian, 
we start with a set of Goldstone fields determined by 
the symmetries of the fiuid. Then we write the most 
general effective action compatible with the symmetries, 
and make use of the conventional derivative expansion 
employed in effective field theories. The main non-trivial 
step in this program is identifying the relevant symme- 
tries; once it is done, the rest is automatic. However, 
the translation from field theory to the conventional lan- 
guage of hydrodynamics requires more work, and can be- 
come quite laborious as one goes to higher orders in the 
derivative expansion. To a large extent, this is related to 
the aforementioned subtleties in the traditional hydrody- 
namical derivative expansion at higher orders. 

Our viewpoint here is that for many purposes it is 
convenient and instructive to consider effective field the- 
ory as a definition of what dissipationless hydrodynamics 
is. As illustrated in the main text, it is straightforward 
then to calculate various physical effects, such as higher 
derivative corrections to the sound wave dispersion re- 
lation. To find what a particular field theory operator 
corresponds to in the traditional language of constitu- 
tive relations may be more challenging (technically, not 
conceptually), but also not really necessary. 

The organization of this paper is as follows. In sec- 
tion [n] we identify the symmetries relevant for describ- 
ing a general perfect fluid at non-zero temperature and 
chemical potential. Here we build on the previous results 
of Refs. H] and |l[, where the effective field theory for- 



malism for zero temperature superfluids and barotropic 
normal fluids was developed. In section IIIII we estab- 
lish the dictionary between field theory and conventional 
thermodynamical variables at the leading order in the 
derivative expansion. 

We continue in section ITVl bv extending the dictionary 
between field theory and hydrodynamics to higher orders 
in the derivative expansion, following the standard pro- 
cedure (see, e.g., Ref. Q). We exemplify this prescription 
by working out explicitly several simple examples. First, 
we consider one-derivative corrections on the field theory 
side. From hydrodynamics, one does not expect any non- 
dissipative terms at this order. In agreement with this 
expectation, we show that on the field theory side these 
corrections can be shifted to higher orders by a field re- 
definition. As a cross-check of our prescription, in the 
Appendix we show that the same procedure of field re- 
definition can be performed in the hydrodynamic theory. 
(In section |V] we prove that the same property holds at 
all orders in the derivative expansion: if a Lagrangian 
term can be removed via a field redefinition, it has no 
effects on the hydrodynamics of the system.) We also 
calculate the correction to the sound wave dispersion re- 
lation following from a sample nontrivial two-derivative 
term. 

We conclude in section IVll bv mentioning a number of 
possible applications of our formalism, some of which we 
are already investigating in detail. 

Our main emphasis here is not on the reformulation 
of the hydrodynamical equations as an action principle, 
but rather on the underlying principles behind this con- 
struction: symmetry and symmetry breaking, Goldstone 
bosons, and derivative expansion. Apart from Refs. [llH] 
mentioned above, earlier works similar in spirit to the 
present paper include (j-Q- Our approach may be use- 
ful for "holographic fluids" where — as emphasized in 
Ref. [ll[ — at low energies the Goldstone dynamics can be 
parameterized without any reference to the microscopic 
theory. 

There has been recent interest in the hydrodynam- 
ics of systems carrying anomalous charges, starting with 
Ref. y. Here, we restrict ourselves to fluids carrying 
ordinary conserved charges, and we avoid the subtleties 
associated with the presence of quantum anomalies. We 
are devoting a companion paper precisely to those sub- 
tleties, and to the resulting interesting effects [12!|. 

Finally, we will deal directly with relativistic hydrody- 
namics, i.e., in our field theory we will impose (sponta- 
neously broken) Lorentz rather than Galilei invariance. 
This choice makes the treatment somewhat simpler. The 
non-relativistic limit can be taken at any stage in our 
analysis, if needed. 



II. FLUIDS WITH CONSERVED CHARGES 

Consider a perfect fluid in d spatial dimensions. Its 
low-energy degrees of freedom can be chosen to be d 



scalar fields 



q^'^^\x,t) 



1=1, 



,d , 



(1) 



giving the comoving (Lagrangian) coordinates of the vol- 
ume element occupying physical (Eulerian) position x 
at time t. This description is reviewed extensivel y in 
Refs. [il 0|j to which we refer the reader for details U8| . 
There is an inherent arbitrariness in labeling the volume 
elements via comoving coordinates. It can be fixed, for 
instance, by choosing these to be aligned with the physi- 
cal ones when the fluid is in equilibrium at some reference 
external pressure. 



(equilibrium) 



(2) 



With this choice of field variables, the fluid's dynamics 
must enjoy the internal symmetries ll| 



a 



a = const 
(f>^ -^R'j4>\ Re SO{d) 
4>' ^ C'(0) , det {de/d<l>-') 



(3) 
(4) 
(5) 



on top of {d + l)-dimensional Poincare invariance. In 
particular, eq. ([5]) corresponds to the fluid's insensitivity 
to (static) non-compressional deformations. Note that 
even though all small perturbations about ([2]) are al- 
lowed, non-perturbatively every field configuration must 
define a time-dependent diffeomorphism between physi- 
cal and comoving space. That is, at any given time it 
must be an invertible function of x |19| . 

Suppose now that the fluid carries a conserved charge. 
There should be an associated U{1) symmetry in our 
field theory. It does not seem sensible to realize this 
symmetry using the (j) fields only: they represent the 
comoving coordinates and physically cannot transform 
under a particle number symmetry. Besides, the (/i^'s are 
non-compact. We should augment the field content to 
represent the particle number symmetry. The most eco- 
nomical addition is a real phase tp{x,t) that shifts under 
it: 



[/(I): 7/)^V + c 



(6) 



We are thus led to construct the low-energy effective 
field theory for the 0^'s and tp, subject to the symmetries 
([S])-® and to Poincare invariance. However there should 
be an additional constraint on the theory, as we know 
that in ordinary perfect-fluid hydrodynamics, the particle 
number current is comoving with the fluid: 

j'" = nuf" . (7) 

(The fluid velocity u^ can be defined through, e.g., the 
entropy current of the fluid.) This guarantees that, in 
normal fluids, sound waves are the only propagating wave 
solutions: the new charge degree of freedom does not in- 
troduce new waves (in contrast to superfluids where there 
are first and second sounds). From our discussion so far. 



it is not obvious how this is going to arise. It turns out 
that enforcing eq. (I7|) is equivalent to imposing a new 
symmetry. If charge flows with the fluid, charge conser- 
vation is obeyed separately within each comoving fluid 
element (recall that diffusion is a dissipative effect and is 
outside the scope of our theory). This means that charge 
conservation is not affected by an arbitrary comoving 
position-dependent redefinition of the charge units. In 
other words, the t/(l) symmetry ^ can be made co- 
moving position-dependent: 



^^^ + /(0/) 



(8) 



where / is a generic function. We will see below that 
this is exactly what we need to enforce ([7]) in our field 
theory. Since the chemical potential will turn out to be 
the simplest invariant under this symmetry, and for lack 
of a better term, we dub this new symmetry 'chemical 
shift'. 

Note that beyond the leading order in the derivative 
expansion, the particle number current is not necessarily 
parallel to the fluid flow. Nevertheless, from the field the- 
ory viewpoint, it is natural to impose the symmetry (|8]) 
to all orders. This automatically implies the existence 
of a quantity conserved along the fiow. This property is 
not obvious in the conventional hydrodynamic language, 
but also the very notion of dissipationless fluid may be 
ambiguous there beyond the leading order. Our view- 
point is that the effective field theory characterized by the 
symmetries ©-(O and ^ provides a natural and un- 
ambiguous definition of dissipationless (nonanomalous) 
hydrodynamics to all orders. 

We are thus looking for the most general relativistic 
Lagrangian that is invariant under ([3])- ([6]) and under ([8]). 
At low energies it should be organized as a derivative ex- 
pansion. At lowest order, we should have one derivative 
per field, because of ^ and ©: 



C^C{d4)',d^l}) 



(9) 



Because of ffl and (O, the 90^ 's should enter in the 
combination [l], 0] 



I 



(10) 



d\ 



^^'-i-'-.^ii-i.0^^^h_,_d^ji._ (11) 



(we define the {d + l)-dimensional e tensor by gOi-'^ = 
-f 1.) J'^ has an important property: it is a vector field 
along which the comoving coordinates do not change: 



J'^d.cl)' = 



1, 



(12) 



Thus it is natural to define the fluid's four-velocity as a 
unit vector aligned with J'': 



1 



-J" 



6=y^Jp^. 



(13) 



We use the 'mostly plus' signature. From a geometric 
viewpoint, J^ is the current of fluid points [l|. 



Now we make use of the chemical shift ([8]) . We notice 
that, because of p^ . the combination 

J'^d^i; (14) 

is invariant. In fact, it is the only invariant at this or- 
der in derivatives. The reason is that the only vector 
that is orthogonal to dnf for generic f{(f>^), is the 'vector 
product' of all the dcj) s, eq. (|lll) . 

To summarize, at lowest order in the derivative ex- 
pansion, the Lagrangian can depend on the 0^'s and on 
tp through J'' and J^d^ij] only. It must be a Poincare 
scalar. We choose to parametrize it as 



S= d^xF{b,y), 



where b is given above, y is defined as 



1 



y = u^^d^ij = - J'^a^V , 



(15) 



(16) 



and F is a generic function. We will see below that F is 
related to the equation of state of our fluid. 

The same field theoretical description of conserved 
charges in hydrodynamics — with emphasis on the same 
symmetry ^ — has been worked out independently by 
Sibiryakov 15]. 

The Noether current associated with t/i's shift symme- 
try © is 



OF 

— u^ = i^,u^ 

dy 



(17) 



and is indeed comoving with the fiuid, as desired. The 
chemical shift ([8|) is an infinite-dimensional symmetry, 
and as a consequence there are infinitely many currents 
associated with it. They are 

J^f) = Fyf{cl>')u'^^f {<!>') r. (18) 

Their conservation is implied by that of j^: 

S^jI^) - f{cb')d^f + Fy djicj,')u'^ ; (19) 



the second term vanishes identically, thanks to p2|) . 

To develop more intuition on why (J15p is the correct 
description of a fluid carrying a conserved charge, let us 
consider how the system behaves in the presence of an 
external gauge field A^. In the field theory, the natural 
way to describe this is to gauge the 'ijj shift symmetry, 
i.e. to replace dfj.tp — >■ df^^' + ^t^ everywhere in the action 
([15]). The resulting action has a non- linear dependence 
on the gauge field. This may appear puzzling from the 
hydrodynamical point of view, where on physical grounds 
one may expect the following linear coupling between the 
gauge field and the fluid 



Oint, — 



fd^xN{(l)')A^J'' = fd^cj)dTN{(l)')u''A^ (20) 



where N{(j)^) is the comoving charge density, and at the 
last step we switched to the comoving frame, i.e. chose 
(p^ as space coordinates and used a?'(0, r) as the dynami- 
cal variables, with r denoting the proper time along the 
comoving worldline. This switch is performed by making 
use of the identity [ij 



dUd-T^ 



dct 



dx^ 



d^x 



1-x^dt^bd^ 



(21) 



To see the relation between the two descriptions let us 
rewrite also the (gauged) fluid action ([TS]) in comoving 
coordinates, 



S= d^(|)dTb-^F{b,^r^J + u^'A^) 



(22) 



The action p2]) depends on -0 only through its time 
derivative, so that the canonically conjugate momentum 



n„ 



b-'K, 



is time independent on the classical solutions, 11^ = 
n^((^^). By making the Legendre transform w.r.t. drip 
(i.e., switching to the Routhian description [20|) we ar- 
rive at the classically equivalent action, 

fd^q)dT [b-^F{b, y{b, n^)) - U^y{b, tt^) + H^m^A^] , 

(23) 
with a linear dependence on Af^. So, the two de- 
scriptions are indeed equivalent upon the identification 
N{(j)) = n^ ((/)), which is consistent with the field the- 
ory/thermodynamics dictionary we are now going to es- 
tablish. 



III. THERMODYNAMICS 

We now want to make contact with the standard hy- 
drodynamic and thermodynamic description of a fluid 
carrying a conserved charge. From the action 



d'^+'^xF{b,y) 



(24) 



we can derive the stress energy tensor by varying with 
respect to the metric. We get 

T^. = [FyV - Fbb)BYJd^^'d,(j)-' + (F - Fyy)7^f,, , (25) 

where Fi, is the 6-derivative of F, and the matrix B^"^ is 
defined as 

B^'^ = a^(/)-fa^0'^ . (26) 

From the definition of J^ it is straightforward to see that 

det B^-^ = b^ (27) 

— which we used to derive T^i, — and 



BYJd^cjj'd,^ 



VfJ-i^ ~r ^/i^i/ — ^/j. 



(28) 



We can thus rewrite the stress-energy tensor in a more 
familiar form: 

T^. = {FyV - Fbb) u^u, -f (F - Fbb)7j^, (29) 

The fluid's energy density and pressure therefore are: 



P = Fylj - F 



p = F- Fbb 



(30) 



Likewise, by comparing pT)) and (O we get that the 
fluid's charge density is 



n = F, 



V ■ 



(31) 



We can get the chemical potential /i, the entropy den- 
sity s, and the temperature T by imposing the thermo- 
dynamics identities 



p + p = Ts + fin , dp — T ds + fi dn 



(32) 



Before doing so, it is worth pointing out that our vector 
J^ is an identically conserved current: 



5m J^ 







(identity) 



(33) 



This follows straightforwardly from its definition. More- 
over, as we already mentioned, it is aligned with the 
fluid's four velocity, J'^ =bu^. These two properties in- 
vite identifying J^ with the entropy current, and b with 
the entropy density, 

s = b . (34) 

That entropy is conserved identically — i.e., 'off-shell' — 
in our non-dissipative field theory, makes perfect sense. 
We could imagine coupling our field theory Lagrangian 
to external sources. These sources could perform work 
on the system, but would not exchange heat with it. In 
such an instance our fiuid would be 'off-shell', but entropy 
would still be conserved |21||. The interpretation of b as 
entropy density has been derived independently in [15|. 
From the first identity in eq. (j31]) , we thus get 



T 



-Fh 



fJ' = y , 



(35) 



which is consistent with the second identity too. (For the 
differential identity, one should express dp in terms of db 
and dFy — dn.) There is of course an ambiguity in the 
overall normalization of s and T — we could multiply s 
and divide T by the same constant, without affecting the 
thermodynamical identities. This is of course related to 
Boltzmann's constant, which does nothing but defining 
the units of temperature. 

In conclusion, our Lagrangian F p3|) is naturally a 
function of the entropy density b and of the chemical po- 
tential y. It can be thought of as a somewhat unusual 
thermodynamic potential: dF — ~Tds + n dfj,. It is re- 
lated to the equation of state p(s, n) or p(T, p,) via either 
of the Legendre transforms in ([3D]) . 



IV. HIGHER DERIVATIVE CORRECTIONS 

Hydrodynamics is naturally organized as a derivative 
expansion. For instance, for the hydrodynamic regime of 
a weakly coupled system of particles, the natural expan- 
sion parameters are the fields' time derivatives times the 
mean free time, and the fields' spacial gradients times 
the mean free path. The standard hydrodynamical and 
thermodynamical variables p, p, u^, etc., correspond in 
our field theory to 'composite operators' involving one 
derivative per (f>^ or ijj field. To reproduce higher-order 
corrections to the perfect fluid hydrodynamics, involving 
derivatives of such variables, we need to include in our 
field theory Lagrangian terms involving correspondingly 
more derivatives than the lowest order ones. Because of 
this, we will number higher derivative corrections start- 
ing from our lowest order Lagrangian. That is, when we 
talk about 'one-derivative terms', we mean terms involv- 
ing overall one more derivative than one per field; and 
so on. Of course, since we will work at the level of the 
action, our field theory will be conservative by construc- 
tion. That is, our approach will not be able to reproduce 
dissipative effects. 

We can adapt to our field-theoretical framework the 
procedure of dealing with higher-derivative hydrodynam- 
ics outlined by Israel and Stewart (IS) Q. The ques- 
tion is — essentially — how to apply thermodynamics to a 
fluid in the presence of spacial gradients and time deriva- 
tives, which typically signal that the system is not in 
complete equilibrium, even though there is some form 
of local thermodynamic equilibrium. IS argue that one 
should proceed as follows. At any given spacetime point 
X, the stress-energy tensor T^^{x) and the charge current 
jfi{x) are perfectly well defined quantities, even for out- 
of-equilibrium systems. For us, they descend straightfor- 
wardly from the Lagrangian, respectively by varying with 
respect to the metric and as the Noether current associ- 
ated with the shift symmetry ^. One defines the local 
energy density p{x) and the local charge density n{x) by 
taking contractions with the local u^{x): 



P 



u'^u-'T., 



-u^J^ 



(36) 



{u'^ is time- like — hence the minus sign.) In the presence 
of gradients, the fiuid velocity field u'^ is itself ambigu- 
ous. For instance, the energy flow and the charge flow 
are typically not aligned with each other. Such an am- 
biguity is harmless 3], and in fact, typically it can be 
used to simplify some algebra. Moreover, p and n are 
particularly well-behaved from this viewpoint, since they 
are unaffected by small variations of u^^ , at first order in 
these variations. One then defines the local values of all 
other thermodynamic variables by applying the equilib- 
rium equation of state to the local p and n thus defined: 



P{x) = Pa{p, n) , p.{x) = p.o{p, n) , 



etc. 



(37) 



The subscript zeroes are there no remind us that we 
should use precisely the same functions of p and n as 



for the fluid in equilibrium, i.e., in the absence of gradi- 
ents. Finally, one goes back to the stress-energy tensor 
and the current, subtracts the perfect fluid part accord- 
ing to the above identifications, and interprets whatever 
is left as the higher-derivative corrections: 



Tf^„ ^ {p + p)ui_,u„ + p rii_,^ + 5T^^ 



Jp = 



Sj^, 



(38) 
(39) 



The IS prescription is in a sense merely a convenient 
definition of what we might want to mean by thermody- 
namical quantities for an out-of-equilibrium fluid. It has 
the advantage of establishing an unambiguous dictionary. 
Moreover, at first order in gradients all such quantities 
are independent of the precise choice of u^^ , thus appear- 
ing perfectly well-defined, and physical. On the other 
hand, at second order and up they all become inherently 
M^'-dependent, and attaching any precise physical mean- 
ing to them becomes more and more dubious [3|]. As an 
important example for us, the entropy density defined as 
above, 



so(p, n) 



(40) 



will not coincide in general with the "observed" entropy 
density — w^s^, where s^ is the entropy current 0]. This 
is well defined even for (slightly) out-of-equilibrium sys- 
tems, because the second law has to hold for them: 
^^s^' > 0. 

As we already emphasized, for us there is no entropy 
production, even off-shell — because our system is conser- 
vative, by construction — and the entropy current is nat- 
urally identified with the identically conserved current 



J'' . 



(41) 



J^ also defines an unambiguous rest frame for the fluid, 
which differs in general from those associated with the 
C/(l) current j'^ and with the stress-energy tensor. We 
will refer to this frame as the 'fleld-theory frame' or the 
'entropy frame'. 

We want to stress that our fleld theory can be taken 
as an independent definition of a non-dissipative fluid 
with mild gradients, that is, mildly out of equilib- 
rium. Its thermodynamics may be ambiguous — a la Is- 
rael and Stewart — but its dynamics are not. For non- 
thermodynamical questions, thermodynamics and the IS 
procedure can be bypassed entirely and the relevant ob- 
servables can be computed directly from the fleld theo- 
retical description. For instance, this is the case for the 
higher-derivative corrections to the sound-wave disper- 
sion law of sect. IIVCI 

We start by considering one-derivative terms in our 
fleld theory. As is well known, absent anomalies [y], 
the only one-derivative corrections one can write down 
for hydrodynamics are dissipative — they are associated 
with shear viscosity, bulk viscosity, and conductivity [lOf . 



As such, they cannot be reproduced by our field the- 
ory. We thus expect that — once interpreted correctly — 
one-derivative corrections to our Lagrangian will be triv- 
ial. As a warmup, we consider one-derivative corrections 
to the dynamics of a zero-temperature superfluid, which 
also turn out to be trivial. 



A. First-order superfluid dynamics 

Consider a relativistic superfluid at zero temperature. 
Its zeroth-order field theoretical description has been 
worked out in [2|. It involves a scalar field ip with a 
shift symmetry "0 — > f/; + a, with a time-dependent vev 
{ip{x)) oc t. The lowest order Lagrangian is 

Co = P{X) , X = {d^f . (42) 

The associated stress-energy tensor and current are 

TO, = -2P'(A)9^Va.V^ + PTit.. (43) 

jO = 2P'{X)d^4, . (44) 

At this order it is natural to define w^ along 9^'(/'- 
„,o _ ^mV- 



(45) 



The minus sign upfront is to some extent a matter of con- 
vention. We are assuming that the superfiuid's ground 
state has ip ex -ft. Taking the appropriate contractions 
with Tq'^ and jg we get 

po = 2P'X -P , no = -2P'yf^X , 



i , Po = P. 
(46) 
These are consistent with the zcro-tcmperature thermo- 
dynamic identities 



p + p = ^j.n 
with chemical potential 



dp = p,dn 



(47) 



(48) 



The function P{X) is thus naturally interpreted as the 
equation of state, giving the pressure as a function of the 
chemical potential (21 : 



Po = P{-yl) 



(49) 



We now add all possible one-derivative corrections con- 
sistent with the symmetries. They take the form 



A/: = G{X)^^'^p^^X 



(50) 



There is another possible structure at this order — 
H{X)\I\ip — which however can be rewritten as above 
upon integrating by parts. This just redefines G(X), 
which is arbitrary anyway. Before proceeding with the IS 
prescription, we notice that A£ can be removed by a field 



redefinition, at the price of introducing higher derivative 
corrections, with two or more derivatives. The reason is 
that it vanishes on the zeroth-order equations of motion, 
for any G{X). Indeed: 



= 9^G(A)F'(A)a^V, 



(51) 
(52) 



where G = J G/P'dX. If we integrate by parts, we get 
precisely the equations of motion associated with the ze- 
roth order Lagrangian (|^^ . As a result A£ can be re- 
moved by redefining tl): 



i,^i,' - \G[X') 



(53) 



As usual, as a byproduct of this redefinition we get higher 
order terms, with two or more derivatives. The effects 
associated with A£ can thus be deferred to higher orders 
in the derivative expansion. Whenever something like 
this happens, the corresponding Lagrangian term is said 
to be 'redundant'. Not surprisingly, if one applies the 
IS procedure to a redundant term before performing the 
field redefinition that removes it, one also gets trivial 
effects at the order under consideration. We prove this 
in sect. |V]in broad generality, and in the Appendix for 
the case under consideration. 



B. First-order fluid dynamics 

For our fluid, the most general first order terms con- 
sistent with our symmetries are 

A£ = /i(6,y) J''d^b + f2{b,y) Jf^d^y , (54) 

where /i and /2 are generic functions. There is in prin- 
ciple another one-derivative structure. 



f3ib,y)rj'^d^J,, 



(55) 



but this is in fact of the same form as the first term in 
([54]) . as J^d^Ju = -\d^b'^. We now show that A£ is 
redundant, and it can thus be removed via a field redef- 
inition. To this end, it is useful to inspect the zeroth 
order equation of motion for ip. It is the conservation of 



d^iFyU^") ^ 



(56) 



Given that J^ is identically conserved, it can be pulled 
out of the derivative. One gets 



J'^df.N = 



N ^ Fy/b 



(57) 



N is the inverse 'entropy per particle', and the above 
equation states the well-known fact that — at zeroth 
order — such a quantity is conserved along the flow. N 
is of course a function of b and y. It is useful to change 



variables in ([54)1 and express everything in terms of b and 

N: 



AC - Mb, N) rd,,b + h{b, N) J^'d^N 



(58) 



The second term vanishes on the zeroth order equation 
of motion ([57l) . The first does not, but we can get rid 
of it via the following trick. We define a new function 
g{b, N) such that 

h{b,N)^dbg{b,N), (59) 

that is, g{b, N) = J /i(6, N) db. We thus get 

AC = J^d^gib, N) + /3(6, N) Jf^d^N , (60) 

where /a = /2 — d^g- The first term now is a total 
derivative — because J^ is identically conserved — while 
the second vanishes on the zeroth-order field equations. 
As a result, all physical effects associated with AC can 
be moved via a field redefinition to higher orders in the 
derivative expansion. As we mentioned, this is consis- 
tent with the absence of non-dissipative one-derivative 
corrections to hydrodynamics. 



C. Sample second order correction 

We leave for future work a thorough study of second 
order corrections in our field theory. Here instead, as an 
illustration, we apply the IS procedure outlined above to 
a sample two-derivative term allowed by all symmetries: 



AC^adf.bd'^b, 



(61) 



where a is an arbitrary coupling constant. Also, for sim- 
plicity we assume that our fluid does not carry any con- 
served charge, in which case the lowest order Lagrangian 
does not depend on our charge field V'l 



Co - Fib) 



(62) 



The contribution to the stress-energy tensor associated 
with AC is 



Ar„^ = -2a dab d^b + 2a bUbP, 



flV 



(63) 

where Pf^i, is the transverse projector of eq. (|28)) . As u^, 
we will still use our zeroth-order definition (IT5|) — i.e. that 
associated with the entropy current J^. The correction 
to the energy density coming from AC therefore is 



Ap = u^u" AT^^ = -2a{u ■ dbf 



(64) 



We are now supposed to derive the other thermodynamic 
variables via the zeroth order relations between p, s, and 
p, which define the equilibrium equation of state. We 
find it convenient to express everything in terms of the 
entropy density. From eqs. (|M| and (PO)) . and setting 
i^y — )■ we get 



p = po - F'{so)As 



(65) 
(66) 



where sq = b is the zeroth-order entropy density. A com- 
parison with (|64p gives immediately 



I 



As = —2a{u ■ dby 



and therefore 



F 



bF" 
Ap = -2a{u ■ dbf 



(67) 



(68) 



Notice that the combination bF" / F' is precisely the 
squared speed of sound [H, Q, so that Ap = c^Ap, as 
implied by our using the zeroth-order equation of state. 
We can now rewrite the full stress energy tensor in 
terms of the 'physical' quantities p, p, s redefined as 
above. We get 



T = T" + AT 



- {p + p)Uf^Ui, +prifj_„ 
2aPf,,,{sDs + 2cl{u ■ dsf) 
aP^.o.P.iid'^sd^s 
2a[u ■ ds) d°'s Pa(fiUy-^ . 



(69) 



The gap in simplicity and clarity between our field the- 
oretical starting point — eq. (|5T|) — and the more standard 
hydrodynamical parametrization of the same second- 
order correction — eq. (15^ — is manifest. For instance, if 
we expand eq. (|6ip in small perturbations about a homo- 
geneous and static background. 



1)^ ^X^ + TT^ 



(70) 



we get directly a correction to the quadratic Lagrangian 
for the sound waves 111 



Co + AC^^ip + p)[7f^~c^,{V-7f)^] (71) 

+ aid^V ■ nf . (72) 

At low momenta, this corresponds to a quartic correction 
to the dispersion law: 



UJ 



ct k^ -f 2a 



(p+p) 



(73) 



Finally, notice what we anticipated above: that the 
entropy density defined following the IS prescription does 
not coincide, at second order, with that associated with 
the entropy current: 



b + As^ -u^J^" 



REDUNDANT COUPLINGS 



(74) 



As we saw in the last section, at the level of our field 
theory certain higher-derivative corrections will be re- 
movable via field redefinitions. We now want to show 
that when this happens, the corresponding corrections to 
hydrodynamical quantities like the current and the stress 



tensor are also trivial. That is, if one is not alert enough 
to detect the possible field redefinitions directly at the 
level of the Lagrangian and goes through the somewhat 
laborious IS procedure, at the end one is left with van- 
ishing corrections to the hydrodynamics of the system 

To see this, suppose that the IS matching has been 
carried out up to n-th order in the derivative expansion. 
One then adds to the action an n + 1-st order term: 



5„+i[$] = S'„[$]+A5„+i[$] 



(75) 



By $ we are collectively denoting all our fields (/>-', tp, 
and in fact the argument we are going to give applies to 
more general situations, e.g. for hydrodynamical systems 
involving more fields. If ASn+i is redundant, there is a 
field redefinition 



$ = $' + G[$'] 



such that 



S^m + ASn+im = Sni^'] + o(a"+2) 



(76) 



(77) 



At n-t- 1-st order, we can drop the 0(9"+^) higher-order 
correction. Clearly, if one applies the IS procedure di- 
rectly to the r.h.s. of this equation, one recovers the n-th 
order results if expressed in terms of physical quantities 
like p, n, etc. — calling the fields $ or $' makes no dif- 
ference from this viewpoint. How does this relate to ap- 
plying the IS procedure to the l.h.s.? Recall that the 
starting point for the IS prescription is the stress-energy 
tensor and the current, and everything else follows from 
there. In our field theoretical framework, these are given 
by functional derivatives of the action with respect to the 
metric and to df^ip, respectively: 



T — —2 



,56' 



SS 



J' 



Sid^.4') 



(78) 



The functional derivative w.r.t. df^ip should be evaluated 
by treating d^ip as a generic vector field, with arbitrary 
variations that vanish at the boundary. This is totally 
unambiguous because our action does not contain undif- 
ferentiated fields. 

Now, the crucial point is that the field redefinition ((76)) 
will necessarily involve derivatives. The reason is that, 
once plugged into Sn , it is supposed to get rid of a term 
involving more derivatives than those contained in Sn- 
That G contains derivatives has two effects: 

1. It mixes the fields with their derivatives; 

2. It mixes the metric with the fields, since the fields' 
derivatives are implicitly contracted via the metric. 

By 'mixing' here we mean a reshuffling of how the action 
depends on the variables involved. 

As a direct consequence of item 2., the stress- 
energy tensor gets 'contaminated' with the equations of 



motion — the functional derivative w.r.t. to the metric ac- 
quires a piece proportional to the functional derivatives 
w.r.t to the fields: 



T = — 2 



SS 



T = 



-2- 



SS 



Sg)"" 



T -2^ 
'^ S^ 



SG 



(79) 
(80) 



where the star denotes the standard integral convolution. 
Therefore, the two stress-energy tensors agree on-shell. 

That the same happens for the current is less trivial 
to see, but equally true. Roughly speaking, it follows 
from item 1. above — the functional derivative w.r.t. d/^,,^ 
acquires a piece proportional to the functional deriva- 
tives w.r.t to the fields — but of course we cannot treat 
"0 and 9^'0 ^-s independent functions in performing func- 
tional variations, so we have to be more precise. A cru- 
cial fact that helps us is the following. Not only does 
the field redefinition (j76|) involve derivatives — it only in- 
volves derivatives. Meaning: the G functional does not 
contain undifferentiated fields. If it did, it would spoil 
the structure of the n-th order Lagrangian — instead of, 
or on top of affecting terms of order n + 1 and above. 
Recall that in our field theory all fields enter the action 
with at least one derivative acting on them. So, for in- 
stance, if we redefined ip SiS tp' + eV''^ and we plugged 
it into So — J F{b,y), we would get corrections to the 
Lagrangian of the form eF'yb and eF' ip u^d^b. The first 
term redefines F . The second does not belong in our 
power counting scheme, because of the undifferentiated 
■0. Perhaps more to the point: in order for the field re- 
definition (1751) to get rid of A5„+i starting from Sn-, it 
has to obey the same symmetries as 5„ and A5„+i, in 
the sense that $' and $ have to transform in the same 
way under these symmetries. Among these symmetries, 
there is shift invariance for all the fields. Therefore, G[^'] 
must be shift invariant. 

We can now compare the currents we would get in the 
$ and $' representations: 



SS 



J 



-5(5,,^) 
SS 



r + 



SS Sido^G) 



sidc^) s{d^,^') 



(81) 
(82) 



The fact that G only involves derivatives of the fields 
allows us to pull the da out of the last functional deriva- 
tive. We can then integrate it by parts (recall that the 
* denotes a convolution), and finally use the fact that 
for a shift-invariant theory, the equations of motion are 
just (minus) the divergence of the corresponding Noether 
currents. We thus get 



J 



f + 



SS 



SG 



S^ Sidf,^') 



(83) 



Like for the stress tensors, the two currents coincide on- 
shell. In the appendix we will carry out the IS matching 



for the redundant coupling ((50)) . and confirm these gen- 
eral results for that case. 

A careful examination of this argument shows that the 
result which we proved has actually nothing to do with 
the IS procedure. Namely, this proof demonstrates sim- 
ply that the on-shell energy-momentum and particle cur- 
rent do not change under field redefinitions. This res- 
onates well with the well-known result that the S'-matrix 
is invariant under field redefinitions [l6| [23|. 

A final comment is in order. Strictly speaking, for 
given Tfj^i, and j^ one can get different results via the IS 
prescription for different choices of w^. So, the findings 
of this section should be interpreted as saying "there is 
a choice of u^^ for which the IS procedure applied to re- 
dundant couplings gives vanishing corrections" . For dif- 
ferent choices of u'^, one gets nontrivial- looking correc- 
tions, which however just amount to the "corrections" 
one would get by boosting the n-th order expressions for 
Tf,u and j^. 



VI. CONCLUDING REMARKS 
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Appendix A: Israel- Stewart matching for first-order 
superfluid dynamics 



Consider the first-order correction (|50p . Its contribu- 
tions to the stress-energy tensor and to the current are 



V = -2G5(^^ d,)X + 2Gatp d,,i^d,xP 




+ ri^,,^C 


(Al) 


j^ = Gd^X-2GUi:d^i^. 


(A2) 



Notice that the tensor structure simplifies considerably 
if one defines a new velocity field 



Let us conclude by mentioning a number of possible 
future directions and open questions. 

Clearly, the major deficiency of this formalism is that 
in its present form it does not allow to discuss dissipa- 
tive phenomena. It appears to be possible to introduce 
these by allowing couplings between the fluid Goldstones 
and an additional soft sector, akin to what happens in 
holographic fluids, where the near-horizon bulk modes 
are playing the role of such a sector fl!] . We leave this 
important challenge for future work. 

Apart from providing a straightforward and clean 
recipe for organizing the derivative expansion, the effec- 
tive fleld theory description brings in other beneflts. Its 
self-consistency implies constraints that are hard to im- 
pose in the conventional hydrodynamical formalism, such 
as unitarity — in the form of absence of ghosts, for exam- 
ple. This may lead to universal inequalities restricting 
the fluid properties, such as the null energy condition (if, 
or an upper bound on anomaly coefficients [12J. 

It is straightforward to extend this formalism to incor- 
porate a larger number of conserved currents. A more in- 
teresting question could be to explore alternative choices 
of symmetries acting on the fiuid Goldstones and to see 
which ones may lead to interesting systems that can be 
realized in nature (some symmetries leading to interest- 
ing systems that are very unlikely to be realized in nature 
have been already explored in studies of massive gravity 

Finally, given the recent interest in hydrodynamics 
with anomalous charges, an obvious application of our 
methods would be to reproduce the associated effects 
via effective field theory. Like for the chiral Lagrangian, 
anomalies at low energies should be encoded in our Gold- 
stone Lagrangian by a Wess-Zumino term. We initiate 
exploring this in a companion paper |12|. 



L " y=X 2P' 
where A^ is a normalization factor: 

A^^ = 1 + AC/{P'X) . 



(A3) 



(A4) 



At first order in derivatives — or in G — this has the ef- 
fect of aligning the full j^ with the velocity field and, 
simultaneously, of diagonalizing the full T^i,: 

Jp = jI + Aj^ = 



= [2P' - 2Ga^P + A£/X)V^Xu 
= {2P'X - 2GXn^ + 2AC)uf,u^ 



rp __ ^0 



(A6) 



Notice also that this redefinition of the velocity field is 
equivalent to the field redefinition ([55]) . 

The corrections to energy- and charge-density associ- 
ated with our one-derivative term are 

Ap = u^'u" AT^^ = G[dip-dX- 20^^ X] (A7) 

An = -u^ Aj^ = G/n/^X [di! ■ dX - 2U%Ij X] (A8) 

We should now define the other thermodynamic variables 
via the zeroth order relations between p, n, n, and p, 
which define the equilibrium equation of state. We find it 
convenient to express everything in terms of the chemical 
potential. From eqs. gg]), (|48]), and ge]) we get 

P^Po- 2P'fioAp (A9) 

P = Po- 2(P' - 2P"pl)poAp (AlO) 

n = no -2{P' - 2P"pl)Ap (All) 
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A comparison with (jA7[) . (jASp gives immediately 

G[d^-dX-2Uil}X] 



Afi ■ 



and therefore 



Ap^ 



2{P' + 2P"X)^/^X 



P'G[dij-dX-2ail;X] 
P' + 2P"X 



(A12) 



(A13) 



Finally, we should express the full current and stress- 
energy tensor in terms of the corrected physical quantities 
defined as above. For the current we have simply 

/' = nu^ (A14) 

— unmodified w.r.t. the zeroth-order one. For the stress- 



energy tensor: 



Tf^i^ = (p + p)Ufj^Ui' + P Vuv 



+ Pm- 2G 



dX ■ dtjj 



p"x 

2P"X+P' 



XD^j 



p' 



2P"X + P' 



(A15) 



where P^^, stands for the orthogonal projector rifj_,y+u^Uiy. 
The second line is proportional to the zeroth-order equa- 
tions of motion 



d" {P'd^,i)) = P'D^ + P"d^'Xd^'4) , (A16) 



and thus vanishes on-shell, as predicted. 
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Earlier works adopting the same parameterization for the 
fiuid degrees of freedom include, e.g., Refs. [13.[l4|. 
To avoid confusion, let us stress that this condition does 
not imply that our variational problem is constrained. 
Locally the condition of invertibility is simply the in- 
equality Aet{d4>/dx) 7^ and does not lead to any con- 
straints on the variations of 0. 

Given a system with Lagrangian C{q,...), where the el- 
lipsis denotes other fields, the fact that the conjugate 
momentum (p) to q is conserved can be used to integrate 
out q. The dynamics is then described by the effective 
Lagrangian CcS = C — pq, where all q dependence should 
be eliminated using p{q, ...) = constant [g]. In our ex- 
ample, tp plays the role of q. Thus —pq — >■ —U^dTip = 
-U^iy - u- A). 

Of course, the same would be true if the entropy current 
were a Noether current, provided the coupling to exter- 
nal sources preserve the corresponding symmetry. The 
difference is that the entropy current would depend on 
the sources in this case. A priori there is nothing wrong 
with this, and this may lead to an alternative dictionary 
between field theory and hydrodynamics. This ambiguity 
may be related to the "integration constants" of anoma- 
lous hydrodynamics [l^ . 

In the interest of full disclosure, we should mention 
that this happened to us for the examples discussed in 
sects. ITV^UVBI 

The two statements are of course related. For instance, 
if one couples the fiuid to dynamical gravity, the rate of 
graviton emission at leading order in Gn is determined 
by the on-shell fluid energy-momentum. 



